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Monte Carlo simulation, experiment and continuum theory are used to examine the anchoring 
exhibited by a nematic liquid crystal at a patterned substrate comprising a periodic array of rect- 
angles that, respectively, promote vertical and planar alignment. It is shown that the easy axis and 
effective anchoring energy promoted by such surfaces can be readily controlled by adjusting the de- 
sign of the pattern. The calculations reveal rich behavior: for strong anchoring, as exhibited by the 
simulated system, for rectangle ratios > 2 the nematic aligns in the direction of the long edge of the 
rectangles, the azimuthal anchoring coefficient changing with pattern shape. In weak anchoring sce- 
narios, however, including our experimental systems, preferential anchoring is degenerate between 
the two rectangle diagonals. Bistability between diagonally- aligned and edge-aligned arrangement 
is predicted for intermediate combinations of anchoring coefficient and system length-scale. 
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I. INTRODUCTION 



Conventional uniform surface treatments for confined 
nematic liquid crystals (LCs), such as rubbed or pho- 
toaligned polymers, are limited to a narrow specification: 
they typically promote vertical or planar alignment where 
both the preferred orientation and associated anchoring 
energy are difficult to alter. In contrast, topographi- 
cally or chemically patterned surfaces permit essentially 
arbitrary control of the easy axis and anchoring poten- 
tial through appropriate adjustment of the pattern fea- 
tures pQ; suitable patterning techniques include atomic 
force microscope scribing of polymer films [2] [3] , micro- 
contact printing of self-assembled monolayers (SAMs) [4]- 
[6] and photolithography [7HT0]. By imprinting a design 
of appropriate symmetry, it is also possible to pattern a 
surface to promote more than one stable alignment orien- 
tation [2j [TTHT4] thus enabling the fabrication of bistable 
devices [2l HQl[T5Hl8] . Very recently, SAM-patterned sub- 
strates have been used to achieve rapid switching [T9] . 

The purpose of the present paper is to investigate the 
alignment behavior of a nematic confined by patterned 
surfaces decorated with a stretched-chessboard-like array 
of rectangles that alternately promote planar and verti- 
cal alignment. This arrangement is of interest since it 
is intermediate between arrays of stripes and squares, 



patterns which have previously been shown to promote 
qualitatively different anchoring behaviors. 

For striped systems, it was found experimentally by Lee 
and Clark [11] that the polar orientation of the bulk ne- 
matic depends on the relative widths of the vertical- and 
planar-promoting stripes while the azimuthal alignment 
consistently lies along the lengths of the stripes. It was 
proposed that the latter was due to elastic anisotropy, 
i.e. the differing energetic costs of various symmetries of 
bulk orientational deformation. Calculations using con- 
tinuum theory [20j [2T] and Monte Carlo (MC) simulation 
of hard particles [22 affirmed this notion and showed 
that it held across the very different length scales over 
which these two theoretical techniques apply. Further, 
it was shown that by adjusting the relative stripe width 
and cell thickness, the polar anchoring angle could be 
altered continuously from planar to vertical [21], 23 . De- 
spite the apparent simplicity of the striped pattern, in 
some scenarios its phase behavior is further enriched: if 
two striped substrates sandwich a cell that is only a few 
molecular lengths thick, the nematic may form separate 
vertical and planar domains "bridging" the film [23] . Al- 
ternatively, if some of the stripes are sufficiently narrow, 
the pattern is neglected by the structure within the ne- 
matic 
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As expected from symmetry [2j [12] , degenerate alignment 
is observed for systems comprising chessboard-like ar- 
rangements of squares promoting competing alignments. 
Experimental studies on SAM-based square-patterned 
systems have shown that, for feature sizes of ~ 30 mi- 
crons, the favored alignment runs along the diagonals of 
the planar-promoting squares [22 . When the patterning 



length-scale is reduced to that accessible to particle-based 
simulation, however, the preferred orientations run along 
the pairs of opposite edges of the squares [26] . An anchor- 
ing transition between these arrangements is then pre- 
dicted for, e.g., a sufficiently weak polar anchoring condi- 
tion. Unlike the striped systems, however, no "bridging" 
behavior has been observed for square-patterned systems, 
even for very thin films. 

Here, we investigate the way in which the strong az- 
imuthal coupling and variable polar anchoring associated 
with nematics on striped substrates segues into the an- 
choring transition found for square-patterned substrates. 
We do this by studying the polar and azimuthal anchor- 
ing behaviors accessible to substrate patternings based 
on rectangles of adjustable length-scale and aspect ratio. 
This is achieved through a combination of molecular-level 
MC simulation, experiment and continuum theory. The 
paper is organized as follows: the MC simulation method- 
ology is described in section II and associated results are 
presented in section III. Some preliminary experimental 
results are then given in section IV and a continuum cal- 
culation is performed and used to reconcile the preceding 
results, in section V. Conclusions are drawn in section VI. 



II. MONTE CARLO MODEL AND 
SIMULATION DETAILS 



obtained for systems of 864 HGO particles of length to 
breadth ratio n = 3, confined between two rectangle pat- 
terned substrates. The substrates were separated by a 
distance L z = 4/^<jo, where <Jo is the particle diameter, 
periodic boundary conditions being imposed in the x- 
and ^/-directions. 

On each substrate, k s was set to a vertical- aligning value 
(k s < 1.0) for a portion of its area and a planar value 
(k s > 2.0) for the remainder and sharp boundaries were 
imposed between the competing alignment regions. The 
patterns on the top and bottom surfaces have been kept 
in perfect registry with one another, as shown in the 
schematic Fig. [2] 
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To initiate this study, we have performed a series of MC 
simulations of rod-shaped particles confined in slab geom- 
etry between two planar walls. The model used is essen- 
tially that described in [26]. Briefly, inter-particle inter- 
actions have been modelled through the Hard Gaussian 
Overlap (HGO) potential [27] which can be seen as equiv- 
alent to the well known Gay-Berne model [28 , stripped 
of its attractive interaction. The particle-substrate in- 
teractions have been modelled using the hard needle- 
wall potential (HNW) [29] where a hard axial needle 
of length (Jok s is placed at the centre of each particle 
(see Figure IT]). The parameter k s provides a molecular- 




FIG. 1: (Color online) Schematic representation of the geom- 
etry used for the hard needle-wall (HNW) particle-substrate 
interaction [29] . 



FIG. 2: (Color online) Schematic representation of rectan- 
gle patterned systems tiled with vertical-inducing (red/dark) 
and planar-inducing (green/light) substrate regions. The az- 
imuthal angle <j> is zero along the x-axis 



Each system has been initialized at low density and gen- 
tly compressed by decreasing the box dimensions L x and 
L y whilst keeping the substrate separation L z and the 
rectangular ratio R = L x /L y fixed. At each density, run 
lengths of 1 million MC sweeps (where one sweep repre- 
sents one attempted move per particle) were performed, 
averages and profiles being accumulated for the final 500 
000 sweeps. 



Analysis has been performed by dividing stored system 
configurations into 100 equidistant constant-z slices and 
performing averages of relevant observables in each slice. 
This yields profiles of quantities such as number density, 
p*(z), from which structural changes can be assessed. 
Orient ational order profiles have also been calculated, 
particularly 
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level control on the surface anchoring properties. The 



results presented in Sections III A and IIIB have been 



which measures variation across the confined films of 
orientational order measured with respect to the sub- 



strate normal. Here N(z) is the instantaneous occupancy 
of the relevant slice. We have also further subdivided 
the system to assess lateral inhomogeneities induced by 
the patterning. Specifically, we have computed profiles 
corresponding to particles residing between the vertical- 
aligning and planar- aligning substrate regions. 



III. SIMULATION RESULTS 
A. Influence of the surface interaction parameter 

In this section, we assess the influence of microscopic 
rectangular patterning on the structure and anchoring 
of a confined LC film. To this end, we present results 
obtained from full compression sequences of MC simu- 
lations performed on a series of systems with differing 
pairs of surface interaction parameters. For reasons of 
space we only show data for high density (p = 0.4) sys- 
tems, and concentrate on the influence of the substrate 
parameters. 

We initially consider systems for which the rectangular 
ratio R = 3. We first used a combination of strong planar 
alignment of the molecules on the substrate (k s = 3) with 
a strong vertical alignment (k s = 0). Then, we slightly 
weakened the vertical alignment (k s = 0.5) before going 
on to use a weak vertical alignment (k s = 1). The high 
density snapshots corresponding to such systems are rep- 
resented in Fig. [3] These clearly indicate that, for all 
combinations of k s values, these systems exhibited a cen- 
tral ordered monodomain at high density. From these 
snapshots it is apparent that the monodomains were all 
aligned in the x-z-plane. This is consistent with what was 
observed for equivalent stripe-patterned systems [23]. It 
suggests, then, that for an edge ratio of 3, rectangle pat- 
terned substrates act rather like well defined stripes, i.e., 
that it is straightforward to pin the azimuthal angle using 
the shape asymmetry of rectangular substrate patterns. 
The second item of note from these snapshots is that 
the surface patterns do not 'bridge' across the film. As 
with square patterns [26], this is presumably due to the 
different interfacial stabilities that would pertain at the 
resultant twist-like and bend-like domain boundaries. 



To substantiate this assessment, we plot, in Fig. [4j the 
corresponding Q zz profiles to quantify the differences be- 
tween these three systems with strong planar anchoring 
(i.e. k s — 3). From these it is apparent that in the ver- 
tical regions, the Q zz value is slightl y grea ter for strong 
vertical parameterizations (see Fig. 



4(a)). This differ- 



ence i s not apparent, though, in the planar regions (see 
Fig. |4(b)[ ). Because of the small differences in these 
graphs, we have also calculated the average bulk tilt an- 
gles. These are presented in Table [I] and confirm that 
these systems exhibit very similar tilt angles in the differ- 
ently anchored regions, typical variations being only 7° 
within each system. These modest variations correspond 




(c) 



(d) 




(e) 



FIG. 3: (Color online) Snapshots of system with L x =3L y and 
different combinations of planar and vertical k s values, (a) 
/c s =0.5; /c s =3; (b) /c s =0; /c s =3; (c) fc a =l; /c s =3; (d) k s =0; 
/c s =2.5; (e) /c s =0; /c s =2 
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TABLE I: Average tilt angle, by surface region, for different 
rectangle ratios L x /L y and substrate couplings k s . 



to nematic monodomains with small undulations of the 
polar anchoring orientation. As expected, on decreasing 
the strength of the vertical anchoring at the surface, the 
bulk alignment becomes increasingly planar: the tabu- 
lated values for Z decrease with increase in the vertical- 
region k s value. Whilst this variation is very weak, the 
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FIG. 4: (Color online) Influence of k s on the p — 0.4 Q zz pro- 
files for the differently-anchored subsystems with rectangle 
ratio R=3 and strong planar anchoring (a) between vertically- 
aligned substrate regions; (b) between planar- aligned sub- 
strate regions 



presence of the vertical- aligning surface regions is sub- 
stantial; the measured tilt angles are far from the 0° that 
would be seen in the absence of these vertical-aligning 
substrate regions. 

We next consider the influence of the strength of the pla- 
nar anchoring regions in the presence of strong vertical 
anchoring regions. This is achieved by simulating k s = 
regions combined with strong planar (k s = 3), moderate 
planar (k s = 2.5) and weak planar(/c s = 2) regions, re- 
spectively. From the variation of the corresponding Q zz 
profiles (see Figs. [5|, it is apparent that decreasing the 
planar anchoring strength, causes the Q zz profiles to be- 
come more positive throughout the confined film. As 
confirmed by the snapshots (Fig. p| and the correspond- 



ing tilt angle data in Table [T[ this means that the cen- 
tral domains become increasingly aligned perpendicular 
to the surfaces. Again, this makes intuitive sense, since 
the relative influence of the vertical anchoring contribu- 
tion would be expected to grow in such a circumstances. 
The bulk tilt angle is, though, more sensitive here to 
change of the planar-region coupling parameter than it 
was previously to change of the vertical-region coupling 
parameter. 
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FIG. 5: (Color online) Influence of k s on the p — 0.4 Q zz pro- 
files for the differently-anchored subsystems with rectangle ra- 
tio R—3 and strong vertical anchoring (a) between vertically- 
aligned substrate regions; (b) between planar- aligned sub- 
strate regions. 



B. Influence of the rectangle ratio R — L x /L y 

In order to assess the influence of the rectangle ratio 
on these small length-scale systems, an equivalent series 



of simulations has been performed on systems patterned 
with R = L x /L y =2 substrate rectangles. Correspond- 
ing high density snapshots are represented in Fig. [7J We 
note from these that again, in all cases, the induced an- 
choring lies in the x-z plane. This indicates that fixing 
the azimuthal angle to coincide with the long edge of the 
rectangle patterns is achievable even with R ~2. The 
influence of the ratio R on the azimuthal angle will be 
investigated in more detail later in the paper. 
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FIG. 6: (Color online) Influence of k s on the p = 0.4 Q zz pro- 
files for the differently-anchored subsystems with rectangle 
ratio R—2 and strong planar anchoring (a) between vertically- 
aligned substrate regions; (b) between planar- aligned sub- 
strate regions. 



Fig. [6] shows the high density Q zz profiles of R=2 sys- 
tems with strong planar anchoring regions and varying 
vertical anchoring strength. These profiles are more sen- 
sitive to the variation of the vertical surface anchoring 
parameter than the corresponding R = 3 were. Specif- 
ically, whilst the R=2 Q zz profiles are similar for k s =0 



and & s =0.5, significantly lower Q zz values were found 
in the central region of the k s =l film. The correspond- 
ing tilt angle shifts are reported in Table |TJ Fig. [8] shows 
the influence of the planar k s parameter on R=2 sys- 
tems with strong vertically aligned regions. Here, the 
Q zz values are strongly influenced by the planar k s val- 
ues, both regions showing increase in Q zz as the planar 
anchoring component is weakened. This corresponds to 
the particles lying increasingly normal to the surfaces, as 
confirmed by the data given in Tab. [I| 
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FIG. 7: (Color online) Snapshots of system with L 
different combinations of planar and vertical k s values. Strong 
vertical anchoring, (a) /c s =0; /c s =3; (b) /c s =0.5; /c s =3; (c) 
fc a =l; k s =3] (d) /c s =0; /c s =2.5; (e) /c s =0; k s =2. 



Whilst the qualitative behaviors of the R =2 and R =3 
systems simulated here were very similar, some quanti- 
tative differences were determined. From the data re- 
ported in Tab. Ill increase in L x /L y was generally as- 
sociated with a decrease in the tilt angles adopted in 
both substrate regions. The one exception to this was 
the substrate combination k s =l k s =3 identified above as 
displaying strong planar character for R = 2. A further 
observation is that the sensitivity of tilt angle to sub- 
strate conditions was consistently greater for the R = 2 



systems than the R = 3. It is not appropriate to examine 
this tilt variation with shape up to the R = inf limit be- 
cause such systems exhibit domain bridging rather than 
a tilted monodomain [23] . 

Finally, we assess the effect of the rectangle ratio on 
the preferred azimuthal angle. In [26 , we saw that the 
azimuthal angle cannot be effectively controlled using 
square-patterned substrates, degenerate alignment being 
seen between the two edge directions. Here, though, for 
both R values investigated, clear orientational pinning 
was apparent. In order to gauge the strength of this pin- 
ning, we plot, in Fig. [9j azimuthal angle distributions for 
substrate-region particles on the different pattern types. 
This confirms a strong departure from the degenerate 
behaviour seen for square patterns. For R = 2, while a 
small but distinct subset of particles aligned along the 
short rectangle edge, long-edge alignment was dominant. 
This effect was even more marked for R = 3 and, in- 
deed, was essentially as strong as that seen for full stripe- 
patterning. 



IV. EXPERIMENT 

To determine the effect of rectangular substrate pat- 
terning at larger length-scales, a set of experimental cells 
was prepared using the microcontact printing technique 
described in [22. In these, SAMs of -COOH and -CF 3 
terminated alkanethiols were used to promote planar and 
vertical alignment, respectively, in the different pattern 
regions. Unlike the simulated systems, however, only one 
patterned substrate was used in these cells; the second 
substrate was always prepared with a uniform vertical- 
aligning treatment. We do not expect the uniform verti- 
cal surface to significantly change the azimuthal align- 
ment since the MC simulations show that the elastic 
distortion is strongly confined to the surface in these 
systems. Moreover, the uniform vertical treatment for 
symmetry reasons cannot itself favor any particular az- 
imuthal orientation. Consequently, these cells comprised 
both vertical- vertical (VV) and planar- vertical (PV) sub 
regions. The cells were prepared with a nominal thick- 
ness 23/im, filled with the nematic material 9CB in the 
isotropic phase and cooled into the nematic phase while 
being observed with a polarizing microscope. 



Fig. 10 shows a polarizing microscopy image of a typi- 
cal cell, one surface of which was patterned with squares 
of S/im width (i.e. a 16 jam period). For this system, the 
maximum transmission in the planar regions occurred 
when the checkerboard was orientated parallel to the po- 
larizers. From this, we deduce that the director had a 
component aligned diagonally across the planar-aligned 
squares. There are two different degenerate states for 
this, corresponding to the two opposite diagonals, which 
are optically identical. Nonetheless, optical consequences 
of splay and twist in the VV regions reveals where bound- 
aries between the differently oriented PV domains lie: at 
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FIG. 8: (Color online) Influence of k s on the p — 0.4 Q zz pro- 
files for the differently-anchored subsystems with rectangle ra- 
tio R=2 and strong vertical anchoring (a) between vertically- 
aligned substrate regions; (b) between planar-aligned sub- 
strate regions. 



an interface between these two the azimuthal angle <j) ro- 
tates gradually by tt/2 from one side to the other, so 
along the diagonal of the square, <j) averages to be paral- 
lel to one of the polariser directions and hence appears 
dark. The VV regions are splayed and twisted in the bulk 
of the cell due to their requirement to accommodate the 
splayed state in the neighboring PV regions. When a 
VV region is surrounded by four PV regions of the same 
hybrid aligned nematic (HAN) configuration, its director 
points in the same diagonal plane as do those of its neigh- 
bors. However, when it is surrounded by PV regions with 
differing HAN configurations, the VV region is not able 
to match its boundary requirements with a simple HAN 
profile. Instead, a twist component develops, which can 
been seen optically as a darker VV square. As these dark 
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FIG. 9: (Online color) Influence of the pattern anisotropy on 
the azimuthal angle distribution of near-substrate particles 
for systems with k s =0 and k s =3. 
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TABLE II: Observed alignment of 9CB on rectangle patterned 
surfaces. 



VV squares only arise where there are neighbors with dif- 
ferent PV region alignments, they form a border around 
domains of the degenerate HAN states. This figure shows 
as well a schematic illustrating these different states and 
their relationship to the microscopy image. This diag- 
onal alignment behavior in VV regions is predicted by 
continuum model calculations [26 if the polar anchoring 
energy promoted by the surface is rather weak. 

A second set of cells was prepared with rectangular 
patterns of different aspect ratios; otherwise, parame- 
ters and conditions were equivalent to those used for the 
square-patterned cell. Microscopy images are shown in 
Fig. pT[a) and (b) and the measured azimuthal align- 
ment angles are displayed in Table \H\ these were obtained 
as for the square pattern by rotating the sample under 
the microscope so as to maximize the extinction between 
crossed polarisers. 

When the symmetry of square checkerboards was bro- 
ken in this way, a number of alignment changes were 
observed. We still have two degenerate azimuthal states 
seen for the regular squares, but the azimuthal orienta- 
tion in the centre is no longer at 45°. On rotation be- 
tween crossed polarisers, dark states were observed for 
anti-clockwise (negative) and clockwise (positive) rota- 
tions. Defining a positive angle (f) from the long axis of 
the rectangle, dark states occur at rotations of -0°, 90-0° 
for one state (</>) and -(90-0)°, 90+0° for the other state 
(-0). However optical degeneracy also has to be con- 
sidered, as the same configurations could be equally ex- 
plained by azimuthal angles of 90-0° and -(90-0)°. Here 
we use the fact that, when observing with parallel po- 
larisers set horizontally, there is no contrast between the 
PV regions and VV regions if the horizontal components 
of refractive index in the HAN and vertical states are 
equal. This allows us to correctly identify the azimuthal 
angles. Optical compensators or direct imaging with flu- 



orescent confocal polarising microscopy (FCPM) could 
also be used to distinguish between these optically de- 
generate states. 



For an aspect ratio of 1.2 at 40° measured from the ver- 
tical. Dark states can be seen on clockwise rotation to 
40° and anticlockwise rotation to 50° , the latter shown in 
fig. [TT[a)(2). Parallel polarizer images also show the dif- 
ference between the bistable states fTTVa) (3). The extinc- 
tion angle is intriguingly close to arctan(l/1.2) = 39.8°, 
the angle of the director if pointing from the center of 
the rectangle to a corner. Additional states are also ob- 
served, as shown in Figure [TT^a)(3), there are regions 
where the director is parallel to the long axis of the rect- 
angle. These regions have a bright state on rotation to 
45° between crossed polarizers. Patterns with an increas- 
ing ratio up to 2 were printed, the final example of which 
can be seen in figure pjjb) with similar behavior to the 
aspect ratio 1.2 case. The observed alignment angles are 
summarized in Table [TTJ Although the MC simulation 
predicts a different alignment direction to that observed 
experimentally, there is a qualitative agreement between 
the two in that the tendency of the nematic is to align 
with the longer side of the rectangle with increasing as- 
pect ratio. 
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FIG. 10: (Color online) (a)Polarizing microscopy image of 9CB aligned by a COOH/CF3 SAM patterned with a checkerboard 
with squares of width 8 /im. The polarizers and analyzers are crossed and are parallel to the checkerboard; (b) Schematic 
showing the bulk director orientation on a square-patterned substrate. A - PV region, B - VV region with identical neighbors, 
C - VV region with different neighboring PV regions. Note the direction of the LC director in a VV region tries to match that 
of its neighbors 



CONTINUUM MODEL 



field 



To reconcile the apparently contradictory observations 
from our MC simulations and experiments, We now ex- 
amine the rectangle pattern using a continuum theory 
approach following the procedure described in our previ- 
ous paper [26 . Our objective is to determine the director 



n(x) = (cos 6 sin (/>, cos cos (/>, sin 0) (2) 

that minimizes the free energy consisting of the Frank 
energy and the surface energy 



F = i / d 3 x K x (V • n) 2 + K 2 [n • (V x n)] 2 + K 3 |n x (V x ii)\ 2 + f dS #(n, n ) 



(3) 



within two simplifying assumptions: that K\ = K% ^ K 2 
with r = K2/K1] and also that only variations in are 
considered, i.e. that the director is confined everywhere 
to a single plane. The first assumption is approximately 
justified for common nematic mesogens, the second is 
reasonable since, as evidenced by Figs. ITT] a) and b), the 
patterned surface promotes variation in the polar coor- 



dinate only and places no restriction on the azimuthal 
orientation. Here, the coordinate system is scaled by the 
thickness of the liquid crystal layer L z , so that the peri- 
ods of the pattern in the x and y directions, L x and L y , 
are dimensionless quantities. 

Here, since L x ^ L y , we require a solution of the form 
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0(x,y,z) = <9 + Yl Yl "TfT ( inm6 UnmZ + B nm e u ™ z )exv[i27r(nx/L x + my/L y )}. (4) 

yL x L y 
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The Euler-Lagrange equation [2^ 




(a)PM images of 9CB aligned on a COOH/CF3 SAM patterned with a checkerboard of rectangles, a = 20 jum, b = 24 fim. 1) A 
uniform region of interest 2) 50° anticlockwise sample rotation with crossed polarizers 3) 50° anticlockwise sample rotation with 

parallel polarizers aligned horizontally. 




(b)PM images of 9CB aligned on a COOH/CF3 SAM patterned with a checkerboard of rectangles, a = 20 /im and b = 40 /im. 1) 
A uniform region of interest 2) 53° clockwise sample rotation with crossed polarizers 3) 53° clockwise rotation with parallel 

polarizers aligned horizontally. 

FIG. 11: (Color online) Polarizing microscopy images of 9CB on a rectangular patterned surface. 



(r cos 2 4> + sin 2 <j>) ^ + (r sin 2 <f> + cos 2 0) |^ + (1 - r) MW)^ + |^ = 



r 



(5) 



is satisfied if the z/™ are chosen 



L X Ly 



2(r + 1) (Llm 2 + L 2 n 2 ) - 2(r - 1) [2L x L y mn sin(20) + cos(2<^) (L^m 2 - ^gn 2 )] . 



r 



(6) 



The remaining coefficients are obtained from the weak 
anchoring boundary conditions, which for the harmonic 
anchoring potential 



g H (0,0e) 



We 



(e-e e f 



are 



±L, 



86 

dz 



(7) 



(8) 



where e (x, y) is the spatially varying easy axis promoted 
by the pattern and the sign corresponds to the direction 
of the surface normal pointing out of the liquid crystal 
layer at the appropriate boundary and where the dimen- 
sionless parameter associated with polar anchoring Lq is 

L > ^ vfc- (9) 

As a guide to interpreting this number, for a cell thick- 
ness L z = 10/im and typical unpatterned values of 
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W = l(T 5 Jm- 2 and ^ w 10p7V (5CB) the value of 
L e = 0.01. 

By inserting Q into ([8| and performing routine calcu- 
lations, the coefficients 9q, A nm and B nm are obtained 



#o = tt/4 

4 - e Cnm (■}()) 

LoVnm{e 1/nm - 1) + (C™ + 1) 

^ m " L flI / nrn (e"»~ - 1) + (e"»». + 1) ' Uij 



where the c nm are the Fourier coefficients of the easy 
axis profile 9o(x,y) at the z = and z = 1 surfaces 
respectively, 



— U>nm — 



yL x L y 



0.0 0.5 1.0 1.5 

Azimuthal Director Angle <j> 



0.5 1.0 1.5 

Azimuthal Director Angle <j> 



n, m odd 
otherwise 



(12) 



FIG. 12: Free energy per unit area as a function of </>, plotted 
for various values of L y /L x . 



Having determined the solution as above, the free energy 
may be evaluated by inserting p] into (p3| and performing 
necessary integrations. The bulk energy is 



E 



LlLli 



m 



B 2 e +v " 



i smh(^ nm j + AA nrn B nm y n 



x { (1 + r) (L 2 X m 2 + L 2 y n 2 ) + (1 - r) [cos(20) (L 2 ra 2 - L 2 n 2 ) + 2L x L y mn sin(20)] } 
- 5Z o^™ [(^nm 6 "^™ + C e+ " nm ) sinh(^ nm ) - 2A nm £ nm i/ n m] . 



(13) 



r 



The surface energy (for each surface) is 
1 



f s — 7T Li x Liyl 10+ / 7 (^nmT%mj V^-nmT-Pnm ^nm)- 

nm 

(14) 
The free energy per unit area as a function of (j) is 
for different values of Lq and L y /L x 
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plotted in Fig. 
with r = 1/2 and L x = 1. For the strongest anchoring 
depicted, where Lq = 0.05, it is apparent that a very 
slight difference in L x and L y is sufficient to break the 
degeneracy of the configurations aligned along the x and 
y axes respectively; once L y /L x > 1.5 there is no stable 
state aligned along the x axis. 

Rectangle-patterned surfaces are surfaces of adjustable 
azimuthal anchoring energy: they may be thought of as 
promoting an effective azimuthal anchoring potential pQ , 
the strength of which is varied by modest changes in the 
aspect ratio of the rectangles. From the plot in Fig.[l2[a), 



it is apparent that over the range of aspect ratios 1 — 2, 
the effective azimuthal anchoring energy varies by a fac- 
tor of roughly 5. The controllability is, however, con- 
tingent on the ability of the pattern to deform the ne- 
matic as measured by the anchoring strength Lq. Free 
energy profiles as a function of (j) for increasing values of 
Lq, corresponding to weaker anchoring, are displayed in 
Fig. [l2|fo-f). As Lq increases, the anchoring transition 
described in our previous paper [26] occurs (Fig. 12 z)) 
whereby the preferred azimuthal orientation, indicated 
by the position of the minimum, is no longer along the 
length or breadth of the square pattern, but lies roughly 
along the diagonal. For the largest value of Lq = 0.2 
plotted, the depth of the effective azimuthal anchoring 
potential is significantly smaller than for the Lq = 0.05 
case. 
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The variation of the preferred azimuthal orientation 
as a function of aspect ratio also depends strongly on 
the anchoring parameter around the anchoring transi- 
tion as may be seen in Fig. 12'c-e), while both the edge- 
(Fig. [l2|a,b)) and diagonally- aligned states (Fig, p^f)) 
are insensitive to the aspect ratio of the pattern . The 
experiment described in the previous section took place 
in this diagonal regime and so the observed alignment 
angle as a function of aspect ratio can be used to esti- 
mate the anchoring parameter; from Tab. [il] we obtain 
a value of Lq ~0.08. Here, the apparent contradiction 
between the Monte Carlo and experimental results is re- 
solved by the continuum prediction that there exist two 
different anchoring regimes characterized by the anchor- 
ing parameter L#; the MC simulations and experiment 
have illuminated the two regimes. In all models the ten- 
dency is for the director to align with the long axis of the 
rectangles. 



VI. CONCLUSIONS 

In this work, a surface prepared with a chessboard-like 
pattern of alternating vertical and planar rectangles has 
been shown to azimuthally align an adjacent nematic liq- 
uid crystal. Moreover, control of the alignment easy axis 
has been demonstrated in both azimuthal and zenithal 
coordinates achieved by adjusting the design parameters 
of the pattern, i.e. the lengthscale, aspect ratio and rel- 
ative anchoring strength of the rectangles. 

Monte Carlo simulations of hard particles predict that 
the preferred alignment direction is along the long edges 
of the rectangles and that only modest aspect ratios 
~ 2 — 3 are required to break the bistable alignment pre- 
viously observed in square-patterned systems [26]. Hence, 
rectangle patterns resemble the striped system exten- 
sively studied [2TJ |23] [24] except that the bulk alignment 
is a monodomain and no bridging behavior, where the 
nematic follows the pattern throughout the film, was ob- 
served; there is only weak modulation in the polar angle 
at the film center, in agreement with that predicted by 
continuum theory. Control of the tilt angle over a range 
of ~ 20° —40° was observed by running a sequence of MC 
simulations adjusting the relative anchoring strength of 
the vertical and planar regions. Simulations with increas- 
ing L x /L y had a distribution of particles increasingly 
strongly peaked around the alignment direction, suggest- 
ing that the azimuthal anchoring strength increases as a 
function of aspect ratio. 

Experimental observations of a nematic aligned on rect- 
angle surfaces prepared by chemical patterning reveal a 
different behavior to that seen in the Monte Carlo simu- 
lation: while for increasing aspect ratio the azimuthal 
alignment direction indeed becomes more oriented to- 
ward the long edge of the rectangles, the alignment di- 
rection is along the diagonals and not the edges. 



The disparity is resolved by our continuum model which 
predicts an anchoring transition from diagonal to edge 
alignment if the anchoring parameter Lq < 0.05 . The 
thin film, only a few particles thick, used in the MC simu- 
lation is inside the edge- aligned strong anchoring regime. 
By fitting the observed variation of the azimuthal anchor- 
ing direction as a function of aspect ratio, we infer that 
the experimental anchoring parameter was Lq ~ 0.08 . 
The proximity of this value to the critical value Lq ~ 0.05 
suggests that with a suitable choice of materials, cell 
thicknesses and patterns that the transition out to be 
observable in a future experiment. 

The key parameter in these systems Lq ~ 0.05, intro- 
duced in eqn.Q, is simply the ratio of the bulk splay 
elasticity to the product of the unpatterned polar anchor- 
ing strength and the film thickness. Thus, it is apparent 
that the edge-aligned behaviour observed in the MC sim- 
ulation part of our study was not dictated by their small 
length-scale - in fact, the small d of the simulated systems 
acted to promote the diagonal alignment regime. That 
the simulations remained firmly in the small Lq ~ 0.05 
regime is, then, a consequence of the very high anchoring 
coefficients pertaining to all models described in terms 
of monodisperse anisotropic particles adsorbed at planar 
substrates. 

Previous studies [30] have shown that the fine-tuning of 
surface anchoring, typically achieved by selection of ma- 
terials and preparation conditions, is necessary to op- 
timize the performance of advanced electro-optic de- 
vices. Patterned substrates of the type considered here, 
however, have the potential to offer a more convenient 
and well-controlled route to achieving such optimisation, 
since all of the surface anchoring parameters, i.e. the 
easy axis, polar and azimuthal anchoring energies, can 
be controlled by adjusting the design features of the pat- 
tern. In such systems, the emergent polar and azimuthal 
anchoring parameters can be viewed as resulting from a 
convolution of the basic material parameters and the im- 
posed patterning. The conventional continuum descrip- 
tion, thus, corresponds to an effective integration of a 
potential spectrum of various surface-printed features on 
length-scales ranging from nanometres to microns. 
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